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Using acoustic methods we have measured nonlinear AC conductance in 2D arrays of Ge-in-Si 
quantum dots. The combination of experimental results and modeling of AC conductance of a dense 
lattice of localized states leads us to the conclusion that the main mechanism of AC conduction in 
hopping systems with large localization length is due to the charge transfer within large clusters, 
while the main mechanism behind its non-Ohmic behavior is charge heating by absorbed power. 
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INTRODUCTION 



Nonlinear hopping conduction in various materials and 
devices was extensively studied, see, e.g., [1, 2] and ref- 
erences therein. Despite relatively large amount of ex- 
perimental and theoretical work aimed at nonlinear phe- 
nomena in static (DC) hopping conduction, the under- 
standing in this research area is far from being complete. 
The common knowledge is that nonlinear effects in hop- 
ping conduction are determined by interplay between the 
so-called field effects (field-induced deformation of per- 
colation paths governing the conduction) and heating of 
the charge carriers. The non-Ohmic effects are usually 
expressed in terms of the dimensionless ratio eEC/kT^ 
where e is the electron charge, E is the electrical field 
strength, T is the temperature, k is the Boltzmann con- 
stant, while C is some characteristic length. Theoretical 
considerations [3] based on different models lead to ap- 
parently similar predictions: non-Ohmic effects become 
noticeable at eEC/kT > 1. However, predictions of 
both magnitude C and its temperature dependence sig- 
nificantly differ. Interplay between the field and the heat- 
ing effects in DC non-Ohmic hopping conduction is es- 
sentially dependent on the relation between the localiza- 
tion length ^ and a typical distance between electrons [1]. 
In two-dimensional (2D) systems with large localization 
length the hopping regime is qualitatively different in 
many respects from the conventional hopping in systems 
with a small ^ [j]. In particular, the nonlinear effects in 
the 2D hopping transport in systems with a large ^ are 
caused by electron heating [1]. 

Much less attention has been paid to non-Ohmic ef- 
fects caused by a high-frequency (AC) field. Most stud- 
ies have been done on modifications of DC conductance 
under electromagnetic radiation. They revealed a rather 
complicated picture of nonlinear behaviors (see [2| and 
references therein), which requires more theoretical ef- 



fort for complete understanding. 



High-frequency hopping conductance is conventionally 
analyzed within the framework of the so-called two-site 
models according to which an electron hops between 
states with close energies localized at two different cen- 
ters. These states form pair complexes, which do not 
overlap. Therefore, they do not contribute to the DC 
conduction, but are important for the AC response. Be- 
ing very simple, the two-site model has been extensively 
studied, see for a review [5-7] and references therein. 
As is well known [6], there are two specific contribu- 
tions to the high-frequency absorption. The first con- 
tribution, the so-called resonant, is due to direct ab- 
sorption of microwave quanta accompanied by interlevel 
transitions. The second one, the so-called relaxation, or 
phonon assisted, is due to phonon-assisted transitions, 
which lead to a lag of the levels populations with respect 
to the microwave-induced variation in the interlevel spac- 
ing. The relative importance of the two mechanisms de- 
pends on the frequency cj, the temperature T, as well 
as on sample parameters. The most important of them 
is the relaxation rate 7(T) of symmetric pairs with in- 
terlevel spacing r^kT. At cj < ^J^kT jh the relaxation 
contribution to the real part of AC conductivity g\{uj) 
dominates. In this case the imaginary component of AC 
condutivity a 2 ^ en, and cri(co') ex a; with logarithmic ac- 
curacy. According to the two-site model, the non-Ohmic 
conductance decreases with the field amplitude [7,, 18,] . 



In this work we study nonlinear AC conductance in 
the samples containing dense arrays of Ge-in-Si quantum 
dots (QDs) using probe- free acoustic method [9]. We will 
show that AC conduction of dense QD arrays is similar 
to that in hopping semiconductors with large localization 
length. 



EXPERIMENT 



Samples We study B-doped arrays of Ge-in-Si QDs 
with densities n = (3 — 4) x 10^^ cm~^ and filling fac- 
tors v ^ 2.5 - 2.85 [10]. The samples were grown by 
Stranski-Krastanov molecular beam epitaxy on a (001) 
Si substrate. The QD array layer was located at 40 nm 
from the sample surface, see the sketch in the left panel 
of Fig. [TJ QDs are shaped as pyramids with a height 
of 10-15 A and a square base having the side of 100-150 
A. The samples were (5-doped with B, the density in the 
doped layer being ~ 10^^ cm~^. In these systems two 
lowest states in Ge QDs are occupied by holes, and the 
third level is partly occupied. The 4th state is split from 
the 3rd one by 18-23 meV 



ll|. Therefore, the 4th level 



is expected to be empty in the temperature domain of 
interest. 





FIG. 1: 

(right). 



Sketch of the sample (left) and experimental setup 



We studied AC conductance of two samples, which 
were annealed and were close in conductance. The aim 
of anneahng was increasing of the sample conductance; 
during the annealing the dots spread and the distance 
between them decreases from initial ~15 nm to almost 
overlapping. The samples studied in [12| had larger con- 
ductance and therefore AC hopping conductance was not 
observed there. 

Procedure We used the probeless acoustic method to 
determine the complex ac conductance, a = <Ji — i(T2, 
from attenuation, F, and velocity, V^ of a surface acoustic 
wave (SAW) induced on the surface of a LiNbOa piezo- 
electric crystal by interdigital transducers, see right panel 
of Fig.dl The sample is pressed to the crystal surface by a 
spring. A SAW-induced AC electric field penetrates into 
the sample and interacts with charge carriers (holes). As 
a result, both F and V acquire additional contributions, 
which can be related to complex a |9|. We single out 
these contributions by measuring the changes induced by 
external transverse magnetic field, AT{B) = F(0) — F(5), 
and similarly AV{B). For these particular samples AV 
is not measurable thus indicating that (72 <C ai. There- 
fore, further we will discuss only AF and ai. Since the 
hole absorption vanishes in very high magnetic fields we 
extrapolate experimentally found AT{B) to B^oo and 
in this way find F(0) = F versus temperature and SAW 
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FIG. 2: Left - Magnetic field dependence of AF at frequency 
143 MHz for different temperatures shown at the curves. 
Right -The same dependences replotted as function of B~^ . 
Sample 1. 



frequency at frequencies 30-414 MHz in the temperature 
domain 1.8-13 K. 

Linear regime Shown in the left panel of Fig. [2] are 
magnetic field dependences of AF at / = 143 MHz and 
5 < 8 T for different temperatures. These dependences 
are replotted in the right panel as functions of 5~^, which 
can be represented as straight lines. The attenuation 
coefficient at 5 = 0, F, is determined by crossing of these 
lines with the axis B~^ = 0. For the situation relevant 
to present experiment the general expression [9] relating 
F with (Ji{uj) can be simplified as 

F (dB/cm) = 4MK'^qA{q) 47rt{q) ai {uj)/SsV . (1) 

Here q = u/V is the SAW wave vector, K is the piezo- 
electric coupling constant in the LiNbOa, A{q) and t{q) 
are geometric factors given in [9|, which depend on the 
dielectric constants of the sample (e^), vacuum (^o) and 
substrate (si), as well as on the depth of the QD layer and 
the clearance a between the sample and substrate. The 
above expression is valid provided (72 <C ai <C esV/4:7rt{q) 
that is the case in the present situation. 

Shown in the left panel of Fig. [3] are temperature de- 
pendences of (Ji for the sample 1 at frequencies 30.1 and 
307 MHz found from Eq. ([1]). The clearance a was de- 
termined from Newton's rings. Frequency dependence of 







307 MHz • 






X 


8x10"' 




/• o 


8x10"' 


T=4.2 K 




^-6X10" 




•/ / 
/ / 


^d ^^''°' 


















>=" S 




/ / 


b" 


/ 




4x10' 


yi 




4x10"' 


/ 






^^ Q ^ 30.1 MHz 




9/ 




2x10"' 


^^^,,,--< 


r 


2x10"' 


" • ' 





1.5 2.0 2.5 3.0 3.5 4.0 
T(K) 



1 00 200 300 

f (MHz) 



FIG. 3: Left - Temperature dependences of g\(u) in the sam- 
ple 1 for / = 30.1 and 307 MHz. a = 5 x 10"^ cm. Right - 
Frequency dependence of g\ in the sample 2 at T = 4.2 K. 
a = 4xl0~^ cm. The lines are guides to the eye. 



F for the sample 2 at T = 4.2 K is shown in the right 



panel. 

Nonlinear regime Shown in the left panel of Fig. |4] are 
magnetic field dependences of AF at frequency / = 253 
MHz at T = 1.8 K for different input SAW intensities. 
Similarly to the linear regime, the values F(0), and ai 
were determined by extrapolation to B^oo. The inten- 
sities W were determined as follows. Firstly, the signal 
having passed through the system was compared with 
the input signal from the generator. This comparison 
allowed determining the transformation losses in both 
transducers plus the losses in input and output trans- 
mission lines (which are considered to be equal). The 
intensity was changed by synchronous changing of the 
calibrated attenuators installed after the generator and 
after the output transducer keeping the summary atten- 
uation (in dB) constant. Knowing the transformation 
losses in one transducer and in the input line we find 
the acoustic intensity after the input transducer. Similar 
measurements were done for / = 30, 143 and 253 MHz. 

Since for each frequency the inequality TL <C 1 (where 
F is the attenuation coefficient in cm~^ while L = 4 — 5 
mm is the sample length) is met, one can neglect the 
difference between the input and output intensities while 
calculating the typical amplitude E of the SAW-induced 
electrical field. It is evaluated from the energy conser- 
vation law 2FW = ai{E)E'^/2. This procedure is valid 
since due to the inequality FL <C 1 higher harmonics of 
the signal are weak. The obtained in this way depen- 
dence cri{E) for sample 1, at T = 1.7 K and different 
frequencies is shown on the right panel of Fig. |4j One 
observe very weak frequency dependence. 

The regime of weak non-Ohmic behavior {E < 100 
V/cm) was investigated in more detail for the sample 2 at 
/ = 142 MHz and different temperatures, T = 3, 3.3, 4.2 
and 7 K. The dependences Acr/cri(0) = [(Ti{E)/(Ti{0) - 
1] on dimensionless parameter {eEC/kT)'^ are shown in 
Fig.O The characteristic length jC is chosen as 7.5 x 10~^ 
cm in order to have A(j/(Ji(0) = 1 at (eEC/kT)'^ = 1. 
One observes an approximate data collapse. 
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FIG. 4: Left - Magnetic field dependences of AF at frequency 
/ = 253 MHz at T = 1.8 K for different input SAW intensities. 
Right - The dependence g\(E^ for T = 1.8 K with the line 
guided to the eye. Sample 1. 
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FIG. 5: Dependence of Acr/cri(0) on the dimensionless pa- 
rameter {eEC/kTf. Here /\a{E) = ai{E) - cri(O), while C 
is a characteristic length chosen as 7.5 x 10~^ cm. / = 142 
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3, 3.3, 4.2 and 7 K. Sample 2. The line is guide 



SIMULATIONS 

Since we are not aware of any analytical theory for 
AC hopping conduction in dense QD arrays we have per- 
formed numerical simulations for a model system con- 
sisting of a L X L square lattice of randomly occupied 
localized states [l3|. We put L = 100 having in mind 
that for this size we know that finite size effects are not 
appreciable in DC conductance [14]. Each of the states 
can be either empty or occupied (double occupancy is 
not allowed). The simulations were performed for the 
filling factor of 1/2, so that the number of electrons is 
half the number of sites. Disorder is introduced by as- 
signing a random energy in the range [— /7, U] to each site 
(in our numerics we have used U = e^ /d where d is the 
lattice constant). Each site is also given a compensating 
charge ve so that the overall system is charge neutral. 
The charges interact via the Coulomb interaction. In the 
following the unit of energy is chosen as the Coulomb en- 
ergy of unit charges on nearest neighbor sites, e^ /d. The 
unit of temperature is then e^ /dk. 

To simulate the time evolution we used the dynamic 
Monte Carlo method introduced in Ref. [15] for simulat- 
ing DC transport. The only difference is that we apply 
an AC electric field E = Eocosut and that, following 
Ref. [16|, we use for the transition rate of an electron 
from site i to site j the formula 



7i j = ^0 



-2r.,7^. 



(e-^-/^^,l) 



(2) 



where Sij is the energy of the phonon and r^j is the dis- 
tance between the sites, tq contains material dependent 
and energy dependent factors, which we approximate by 
their average value; we consider it as constant and its 
value, of the order of 10~^^ s, is chosen as our unit of 
time. Consequently, uj is measured in units of Tq^ while 
the electric field is measured in units of e/d'^ . The con- 
ductance (in Siemens) is expressed in units of 47reo(i/ro. 
Note that our model is oversimplified comparing to real- 
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FIG. 6: AC conductivity as function of frequency (left panel) 
and of electrical field (right panel) for T = 0.03. 



istic dense QD arrays. In particular, the model assumes 
a hydrogen-like wave function with localization length 
^ and independent of the phonon wave vector electron- 
phonon couphng, we ignore dielectric susceptibility of the 
host, etc. Therefore, we hope to reproduce the behavior 
of AC conductance only qualitatively. However, we be- 
lieve that the model reproduces main physics of the dense 
arrays - interplay between Coulomb correlations and dis- 
order as well as relatively large ratio between ^ and a 
distance between the dots. To decrease the number of 
independent parameters we put (, = d. 

Shown in Fig. [6] (left panel) is the conductance ver- 
sus frequency for T = 0.03 (corresponding to the Efros- 
Shklovskii regime) and ^o = 0.003 which we know is close 
to the upper limit of the Ohmic regime for DC trans- 
port [l4J. Dependences of a on the field amplitude for 
different frequencies are shown in the right panel. As it 
is seen, the conductance grows with frequency, as well 
as with electric field that qualitatively agrees with the 
experimental results shown in Figs. [3l H] (right panels). 

We can present the results in a different form: assum- 
ing that the non-Ohmic behavior is due to heating of the 
charge carriers we would expect that the conductance at 
a large field will be the same as the Ohmic conductance 
at a higher temperature corresponding to the electron 
temperature. This is known to be approximately true 
for DC conductivity [IJ, ll7|. Thus, we expect a relation 
of the form 



a{Tef^,E,uj) = cr(T,0,a;) 



(3) 



where Teff(£^) is the effective electron temperature while 
cr(T, 0,6^;) is the Ohmic AC conductance. In the sim- 
ulations we can check the accuracy of this relation by 
independent calculation of cr(a;) and Teff. The latter was 
found from direct fitting of the Fermi function to the 
actual distribution function of occupied sites obtained 
by averaging over 20 states. We also found the Ohmic 
conductivity as function of temperature for different fre- 
quencies. If relation (J3]) holds we should find at a given 
frequency the following: as the electric field is increased, 
both the effective temperature and the conductivity in- 
crease in such a way that the points a{T^E^uj) follow 



the same curve as the Ohmic conductivity. We plot a 
versus l/^/T^ for different frequencies in Fig. [7] (left 
panel). As we can see, relation (J3j) is approximately sat- 
isfied. Therefore, we have used it as a definition of the 
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FIG. 7: Left - Conductance as a function of 1/V^eff for var- 
ious frequencies. The large symbols connected by lines show 
the ohmic conductivity (for which Teff = T). The small sym- 
bols show the non-Ohmic conductivity at a finite field. Each 
symbol corresponds to a different frequency. Right - Rela- 
tionship between g(T,E)E^ and T^^ - T^ for sample 1 at 
T= 1.8 K. 



electron temperature allowing us to extract it from ex- 
perimental data on <j(T, £^,cj) and (t(T, 0,6^;). Shown in 
the right panel of Fig. [3 is the relationship between so- 
defined Teff and the absorbed power, (j{T^E^uj)E^^ for 
the lattice temperature of 1.8 K. Here we used averaged 
data for different frequencies since the frequency depen- 
dence is rather weak. The results are compatible with 
the dependence T^^ - T^ oc a{T,E)E'^. Interestingly, 
a similar dependence was reported for amorphous InO^^ 
films, which are on the dielectric side of a insulator-to- 
superconductor transition [18|. Though such a depen- 
dence was predicted for a metallic state [19], we are not 
aware of any analytical theory explaining such a depen- 
dence in the insulating phase. 

Conclusion Using acoustic methods we have mea- 
sured real part of nonlinear AC conductance in insulat- 
ing 2D arrays of Ge-in-Si QDs. The observed features - 
inequality ai ^ cr2, increasing of ai with the field am- 
plitude, its relatively weak frequency dependence, large 
energy exchange length C > I/a/^ ~ essentially contra- 
dict to the predictions of the two-site model discussed in 
SeCn The dissipation is rather due to electron transitions 
in multi-site clusters. This is not surprising because the 
two-site model requires the close pairs to be isolated that 
can be met only in dilute systems. The observed non- 
Ohmic behavior is attributed to heating of the charge 
carriers by AC electrical field, as in the DC case [1[. We 
have also performed simulations of AC conductance in a 
dense 2D lattice of localized states taking into account 
both distribution of their energies and Coulomb corre- 
lations. The results qualitatively agree with experiment 
confirming the above conclusion. 
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